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FINAL  REPORT  EVALUATION 
VLSI  Design  For  Reliability  -  Current  Density 

The  growth  in  the  use  of  custom,  semi-custom,  and  application  specific  ICs  (ASICs)  in  AF 
systems  has  brought  several  changes  to  the  reliability  qualification  process.  This  occurred 
because  the  costs  associated  with  reliability  characterization  has  become  a  major  part  of  the 
device  cost  when  there  are  only  a  few  hundred  or  thousand  of  such  devices  produced.  Further, 
the  cost  of  a  second  or  third  pass  at  the  design/fabricate/test/evaiuate  reliability  cycle  is  very 
time  consuming  as  well  as  expensive. 

Analyzing  an  1C  design  for  susceptibility  to  life  limiting  failure  mechanisms  is  one  approach  to 
enhance  the  probability  that  at  the  end  of  the  first  pass  of  the  design/fabricate/test/evaluate 
reliability  cycle,  the  device  will  have  an  acceptable  reliability. 

This  effort  describes  a  more  accurate  and  general  approach  to  calculating  the  expected  value  and 
variance  of  the  currents  in  a  digital  CMOS  circuit,  enhancements  to  extraction  of  the  metal  RC 
networks  and  circuit  descriptions  appropriate  for  probabilistic  simulation,  and  the  computation 
of  current  density  waveforms  in  the  bus  metal.  These  calculations  address  susceptibility  to 
electromigration,  and  may  be  extendible  to  power  estimation  and  IR  drop  problems. 

Other  efforts  sponsored  by  Rome  Lab  in  the  area  of  design  for  reliability  include  simulation 
methods  for  hot  electron  degradation  in  CMOS  ICs,  where  realistic  charge  distributions  are 
employed  in  modeling  the  localized  oxide  damage. 

'yrtsMOJkr^' 

MARTIN  J.  WALTER 

Design  &  Diagnostics  Branch 

Microelectronics  Reliability  Division 
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DESIGN  FOR  RELIABILITY-CURRENT  DENSITY 

Final  Report 

Contract  F30602-88-D-0028(Task  N-0-5700) 

Abstract 

This  report  describes  the  work  accomplished  during  the  period  April  1,  1990,  to  July  31, 
1991,  on  reliability  analysis  of  VLSI  CMOS  circuits  for  electromigration  prevention.  The 
work  is  a  continuation  of  research  on  the  same  subject  that  was  supported  under  the  same 
contract  during  the  previous  year.  The  main  emphasis  of  the  work  is  the  computation  of  the 
expected  value  and  variance  current  density  waveforms  in  sections  of  the  bus  for  estimating 
mean-time-to-failure  (MTF)  due  to  electromigration  effects,  as  well  as  to  determine  the  width 
of  the  bus  sections  in  order  to  keep  the  MTF  above  a  certain  safe  limit.  The  work  involved 
three  subtasks:  (1)  Further  development  and  enhancement  of  the  probabilistic  simulation 
method  to  compute  the  current  and  voltage  waveform  statistics;  (2)  development  of  hierarchi¬ 
cal  extraction  software  for  extracting  transistor  circuit  description  and  parasitics  from  hierarch¬ 
ically  specified  CIF  layout  specifications;  and  (3)  computation  of  the  expected  value  and  vari¬ 
ance  current  density  waveforms  in  the  RC  model  of  the  bus,  given  the  expected  value  and 
variance  current  waveforms  at  the  bus  contacts. 

In  probabilistic  simulation,  the  statistical  approach  for  computing  the  current  and  voltage 
statistics  has  been  refined  and  made  more  accurate  and  more  general,  especially  in  the  deriva¬ 
tion  of  the  variance  of  the  current  waveform  and  the  expected  value  of  the  gate  delays.  The 
new  derivations  are  being  implemented  in  a  new  hierarchical  probabilistic  simulator. 

In  the  area  of  layout  extraction,  a  hierarchical  extractor  has  been  developed  and  linked  to 
the  OCT/VEM  design  data  base  system.  The  extractor  accepts  transistor  artwork  circuit 
specified  in  CIF  in  hierarchical  format  and  produces  a  circuit  description  with  the  same 
hierarchical  format  and  produces  a  circuit  description  with  the  same  hierarchy. 

In  the  subtask  on  power  and  ground  bus  analysis,  a  numerical  approach  has  been  imple¬ 
mented  to  compute  both  the  expected  value  and  the  variance  current  density  waveforms  in  the 
bus  for  MTF  estimation  of  electromigration  effects.  Previous  methods  have  concentrated  only 
on  the  expected  value  or  considered  only  dc  components  of  the  current  waveforms.  The 
approach  uses  network  reduction  and  sparse  matrix  techniques  for  fast  solution. 
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DESIGN  FOR  RELIABILITY-CURRENT  DENSITY 

Final  Report 

Contract  F30602-88-D-0028(Task  N-0-5700) 

1.  Introduction 

This  task  is  concerned  with  the  development  of  computer-aided  design  techniques,  as 
well  as,  software  tools  for  reliability  analysis  and  design  of  Application  Specific  Very  Large 
Scale  Integrated  (ASIC  VLSI)  CMOS  circuits.  There  are  many  failure  mechanisms  that  have 
to  be  considered  when  analyzing  the  reliability  of  integrated  circuits.  These  include  metal  and 
contact  electromigration,  power  supply  voltage  drop  and  noise,  transistor  degradation  induced 
by  hot-electron  injection  into  the  oxide,  oxide  breakdown,  electrostatic  discharge  through  input 
protection  circuit,  system-noise  induced  latchup,  and  charged-particle-induced  soft  errors. 

Since  wearout  mechanisms  leading  to  poor  reliability  are  a  result  the  circuit  behavior 
over  typically  long  periods  of  time,  they  depend  on  the  combined  effects  of  a  large  number  of 
possible  input  waveforms  rather  than  worst-case  conditions.  One  of  the  main  issues  addressed 
in  this  task  is  the  development  of  efficient  techniques  for  representing  the  combined  effects  of 
all  possible  input  waveforms  by  one  or  possibly  two  waveforms  (expected  and  variance 
waveforms),  and  the  derivation  and  implementation  of  fast  techniques  for  generating  these 
waveforms.  This  has  led  to  the  development  of  a  new  methodology  which  we  call  probabilis¬ 
tic  simulation  [1-2]. 

In  this  task,  we  concentrate  our  efforts  on  a  specific  reliability  problem,  namely,  elec¬ 
tromigration  (EM).  In  [2],  we  showed  that  MTF  due  to  electromigration  depends  on  the 
expected  value,  or  average,  as  well  as  the  variance  current  density.  Probabilistic  simulation 
provides  a  very  efficient  way  for  obtaining  these  current  density  statistics  without  the  need  for 
having  to  carry  out  exhaustive  simulation. 
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In  our  first  annual  report,  we  described  in  some  detail  our  initial  work  on  the  probabilis¬ 
tic  simulation  approach  and  its  application  to  compute  current  density  and  estimate  MTF  due 
to  electromigration.  In  this  second  annual  report,  we  describe  our  accomplishments  during 
this  past  year.  In  section  2,  we  summarize  our  continued  work  on  the  probabilistic  simulation 
approach.  In  section  3,  we  describe  the  work  on  hierarchical  extraction;  namely,  the  extrac¬ 
tion  of  hierarchical  transistor  circuit  specifications,  including  parasitics,  from  hierarchically 
specified  layout  CIF  files.  In  section  4,  we  report  on  the  development  of  sparse  matrix  tech¬ 
niques  for  the  computation  of  density  waveforms  in  the  bus  given  the  expected  value  and 
variance  current  waveforms  at  the  contacts  to  the  bus.  In  section  5,  we  summarize  our  results 
and  discuss  future  work. 
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2.  Probabilistic  Simulation  Techniques 

Simply  stated,  probabilistic  simulation  attempts  to  compute  current  and  voltage  statistics, 
in  the  form  of  average  or  expected  value  and  variance  waveforms  in  a  circuit  for  all  possible 
input  excitations,  without  the  need  of  having  to  perform  exhaustive  simulation.  Last  year,  we 
reported  on  the  derivation  of  this  new  probabilistic  simulation  technique  and  its  implementa¬ 
tion  in  program  CREST  [1].  During  this  year,  we  have  embarked  on  the  development  of  a 
new  hierarchical  simulator  in  order  to  be  able  to  analyze  very  large  designs.  In  contrast, 
CREST  accepts  only  flattened  circuit  descriptions.  The  new  simulator  keeps  a  hierarchical 
data  structure  that  will  allow  hierarchical  probabilistic  simulation  techniques  to  be  imple¬ 
mented.  In  addition,  we  have  proved  and  implemented  new  derivations  for  the  probabilistic 
simulation  technique.  The  new  derivations,  which  are  explained  below,  extend  the  previous 
ones  and  makes  them  more  accurate.  We  first  give  a  brief  review  of  probabilistic  simulation. 
The  details  of  the  new  derivation  and  its  implementation,  together  with  comparisons  with 
CREST  results  and  with  exhaustive  SPICE  simulation  for  a  number  of  typical  examples  are 
included  in  a  M.S.  thesis  report  [3], 

2.1  Probabilistic  Waveforms 

Probabilistic  simulation  can  be  described  as  an  event-driven  simulation  related  to 
switch-level  simulation.  As  described  in  CREST  [1]  one  must  stop  short  of  describing  the 
input  signals  as  stochastic  processes.  Therefore  the  probabilistic  waveform  description  is 
chosen  as  the  description  of  choice,  which  is  somewhere  between  a  deterministic  waveform 
and  its  stochastic  counterpart.  Each  input  is  characterized  by  the  probability  that  it  is  high, 
which  means  that  the  corresponding  node  has  high  logic  value.  Also  the  times  where  switch¬ 
ing  occurs  are  given  as  deterministic  variables.  This  description  is  not  complete  though.  To 
establish  a  complete  description  of  the  input  at  every  time  point,  the  transition  probabilities 
must  be  included  in  the  description.  Actually  only  the  probability  of  one  transition  is  required 
(we  will  use  the  probability  that  the  input  switches  from  low  to  high  Pih  in  the  rest  of  this 
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thesis),  because  given  Ph(t  ),  Ph(t+),  and  Pu,(t),  then  Pw(t)  can  easily  be  calculated  from  the 
relation 

Ph(t+)  ~  Ph(0  =  Pih(t)  -  Pw(t)  (2.1) 

In  this  way  using  probabilistic  waveforms  we  can  describe  a  large  number  of  inputs.  An 
example  of  a  probabilistic  waveform  is  shown  in  the  following  figure.  Using  this  example  the 
use  of  a  probabilistic  waveform  will  be  explained. 


0  1  2  3  4  time 


Figure  2.1.  Example  of  probabilistic  waveform 

In  Figure  2.1  the  numbers  over  the  arrows  give  the  probability  that  a  signal  will  switch 
from  low  to  high,  while  the  numbers  on  the  flat  part  give  the  probability  that  the  signal  is 
high. 

Now  let  us  consider  the  time  point  t  =  1.  The  probability  that  the  input  is  high  before 
t  =  1  is  0.5,  while  the  probability  of  being  high  after  t  =  1  is  0.75.  The  transition  probabil¬ 
ity  P|),(t)  at  that  time  point  is  0.3.  Applying  the  above  formula  in  (2.1)  the  probability  for  a 
transition  from  high  to  low  is  Pw(t)  =  Plh(t)  -  Ph(t+)  +  Ph(0  =  0  05. 
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2.2  Current  Waveform  Description 

As  has  been  stated  above  the  objective  is  to  derive  the  characteristics  of  the  stochastic 
process  describing  the  current  through  the  bus.  One  complete  description  of  the  process 
would  be  to  evaluate  the  current  induced  by  every  input  vector,  along  with  the  probability  of 
this  taking  place.  Although  this  is  possible  (because  we  are  dealing  with  a  finite  space)  it  can 
prove  to  be  very  time  consuming.  Consequently  we  have  to  limit  ourselves  to  more  incom¬ 
plete  descriptions.  One  such  description  involves  the  expected  value  waveform  and  the  vari¬ 
ance  waveform  (or  R(t,t)  of  the  corresponding  stochastic  process). 

The  expected  value  waveform  represents  the  average  value  of  all  the  paths  at  each  time 
point.  One  can  visualize  it  as  the  weighted  sum  of  all  the  currents  drawn  by  the  circuit. 

2 

The  variance  waveform  is  calculated  as  V(i(t))  =  E[(i(t)  -  E[(i(t))])  ]  at  every  time 
point.  As  it  is  becoming  evident  we  will  approximate  the  distribution  of  the  current  at  each 
time  point  with  a  gaussian  distribution.  This  is  a  simplification;  but  the  accuracy  of  the 
results  has  validated  our  assumption. 

The  relationship  between  the  expected  value  and  the  variance  waveforms  and  MTF  has 
been  established  in  [2].  The  remainder  of  this  section  will  be  on  how  to  derive  accurately  and 
efficiently  these  two  very  important  quantities  given  the  input  waveforms  and  the  circuit 
topology. 

2.3  Simulation  Algorithm 

The  input  waveforms  given  in  the  probabilistic  waveform  format  are  considered  given 
and  so  is  the  description  of  the  circuit.  The  objective  is  to  calculate  the  waveform  at  the  out¬ 
put  node  and  the  expected  value  and  variance  waveforms  of  the  current  through  the  gate.  As 
it  is  known  that  transitions  cause  output  events  the  simulation  proceeds  in  an  event  driven 
fashion. 
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The  circuit  is  first  partitioned  into  an  interconnection  of  channel-connected  subcircuits. 
Simply  stated,  given  a  probabilistic  representation  of  the  voltage  waveforms  at  the  gates  of  the 
transistors  in  a  subcircuit,  the  transistor  conductances  in  the  subcircuit  become  random  vari¬ 
ables.  Graph  reduction  techniques  are  then  applied  to  derive  an  equivalent  probabilistic 
macromodel  of  the  subcircuit  with  lumped  output  capacitances.  The  macromodel  is  then  used 
to  compute  the  output  voltage  probability  waveform,  the  expected  value  and  variance  current 
waveforms  drawn  by  the  subcircuit,  as  well  as  the  delay  through  the  subcircuit.  Probabilistic 
simulation  has  been  implemented  in  CREST  [1]  and  its  validity  has  been  established  by  simu¬ 
lating  large  industrial  examples. 

A  key  to  the  accuracy  of  the  results  of  the  probabilistic  analysis  of  a  subcircuit  is  the 
accuracy  of  the  graph  reduction  step  in  the  process  of  obtaining  an  equivalent  macromodel. 
We  derive  an  alternative  reduction  technique  that  is  shown  to  be  superior  to  the  previous  one, 
especially  in  the  computation  of  the  variance.  The  accuracy  of  this  new  method  has  been 
established  by  comparing  the  results  to  those  obtained  through  exhaustive  analysis. 

2.4  Single  Gate  Reduction  Algorithm 

The  analysis  of  a  single  channel  connected  partition  of  a  circuit  proceeds  as  follows.  For 
a  fully  complementary  gate  we  separate  the  p-part  and  the  n-part  and  simulate  them 
separately.  Then  we  map  the  part  we  are  simulating  onto  a  graph  as  shown  in  Fig.  2.2,  thus 
producing  a  simple  non-directed  graph. 
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Figure  2.2  Mapping  onto  a  graph 

The  objective  of  the  whole  process  is  to  reduce  the  graph  on  the  right  to  a  single  edge  as 
shown  in  Fig.  2.3. 


Figure  2.3  Reduction  to  an  equivalent  edge 


In  this  part  the  procedure  through  which  we  derive  the  equivalent  conductance  of  the  cir¬ 
cuit  and  the  probabilities  at  all  the  necessary  time  points  will  be  explained  in  detail. 

We  will  consider  only  fully  complementary  CMOS  circuits;  thus  the  subcircuit  can  be 
further  partitioned  into  a  p-part  and  an  n-part.  Since  these  two  parts  are  complementary,  we 
will  consider  the  analysis  of  the  n-part  only.  Each  transistor  in  the  partition  has  a  conduc- 


10 


tance  associated  with  it  when  it  is  ON  with  value  g^,.  When  the  transistor  is  OFF  it  will  have 
zero  conductance.  The  probability  of  an  nMOS  transistor  being  ON  is  the  probability  of  its 
gate  node  being  "high"  while  the  probability  of  a  pMOS  transistor  being  ON  is  the  probability 
of  its  gate  being  "low."  Thus  the  conductance  of  a  transistor  can  be  viewed  as  a  random  vari¬ 
able.  The  expected  value  and  the  variance  of  the  conductance  are  defined  respectively  at  any 
time  t  other  than  switching  time  as  follows: 


E[g]  =  gon  X  Ph(0 

(2.2) 

V(g)  =  gon  X  Ph(0  ~  (gon  X  Ph(t))2 

(2.3) 

As  stated  before  we  aim  at  finding  an  equivalent  edge  between  the  output  node  and  the 
ground  node  (Vdd  node  for  the  p-part),  so  that  we  can  calculate  the  current  drawn  at  the 
ground  contact.  We  assume  for  the  time  being  that  the  circuit  being  eliminated  can  be  elim¬ 
inated  in  a  series-parallel  fashion,  which  is  the  most  prevalent  case.  The  probability  that  the 
voltage  at  the  output  node  is  high  must  also  be  calculated.  We  will  consider  the  reduction  of 
two  basic  cases:  one  is  when  two  transistors  are  in  parallel  and  the  other  when  two  transis¬ 
tors  are  in  series.  Assuming  independent  inputs,  the  parallel  case  becomes  straightforward. 

We  will  use  the  following  notation  in  the  sequel: 

x,y:  the  conductances  of  the  two  edges  (transistors)  to  be  reduced  to  one 
g:  the  equivalent  conductance 
P^:  probability  that  x  is  conducting 
Plhx:  probability  that  x  is  switching  OFF  to  ON 

It  is  also  evident  that  when  the  equivalent  edge  is  ON  the  output  node  is  "low,"  since  we 
are  considering  the  n-transistor  block. 

In  the  parallel  reduction  case  (shown  in  Fig.  2.4)  the  equivalent  conductance  is  given  by: 


g  =  x  +  y 


(2.4) 
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Figure  2.4  Parallel  reduction 

The  stochastic  quantities  describing  the  equivalent  edge  are  given  by  the  following 
expressions,  where  the  subscript  p  (for  parallel)  denotes  those  quantities  associated  with  the 
equivalent  edge.  These  results  follow  from  simple  probability  theory  and  their  proof  is  omit¬ 
ted. 


^hp  —  ^hx  ^*hy  ^hx^hy 

(2.5) 

^Ihp  =  Pflix^ly  ^Uiy^Ix  —  ^Ihx^Diy 

(2.6) 

Efg]  =  E[x]  +  E[y 1 

(2.7) 

V(g)  =  V(x)  +  V(y) 

(2.8) 

These  equations  are  not  approximations  but  exact  expressions  of  the  quantities  for  the 
equivalent  edge. 

The  series  case  (shown  in  Fig.  2.5)  is  not  as  straightforward. 


Figure  2.5  Series  reduction 

The  probability  that  the  equivalent  edge  is  high  and  the  probability  that  it  switches  are 
propagated  by  the  following  expressions  which  again  follow  from  elementary  probability 
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theory  (the  quantities  associated  with  the  equivalent  edge  are  here  denoted  by  the  subscript  s). 

Phs  =  ^hx^hy  (2.9) 

^lhs  =  Plhx^hy  ^Lhy^hx  ~  Plhx^lhy  (2.10) 

The  above  equations  are  also  exact.  But  when  it  comes  to  propagating  the  expected  value 
and  the  variance  of  the  equivalent  edge  we  have  to  resert  to  some  approximation.  In  [1]  the 
following  expressions  have  been  proposed: 


E[g]  E[x]Phy  +  EfyJP^ 

V(g)  _  V(x)  ,  V(y) 
E[g]4  E[x]4  E[y]4 


(2.11) 

(2.12) 


These  expressions  give  only  a  rough  approximation  of  the  values  in  question.  They  can 
be  greatly  improved.  As  proposed  in  (8]  a  Taylor  expansion  of  the  expression  (2.13)  for  the 
conductance  of  the  equivalent  edge  for  the  deterministic  case  is  obtained  and  by  the  appropri¬ 
ate  algebraic  manipulation  a  better  approximation  can  be  obtained.  So  starting  from  the  fol¬ 
lowing  expression 
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V(g)  = 


E[y] 

E[x]  +  E[y] 


I4 

V(x)  + 


E[x] 

E[x]  +  E[y] 


V(y) 


(2.17) 


Unfortunately  these  expressions  do  not  give  accurate  results  especially  when  Ph  is  rather 
small.  But  they  are  very  important  because  the  calculation  of  the  expected  value  and  the  vari¬ 
ance  is  decoupled  from  the  correct  propagation  of  the  probability  at  each  node.  This  can  be 
helpful  in  cases  where  correlation  hampers  the  propagation  of  the  probabilities.  But  the  fact 
that  it  cannot  be  used  reliably  in  cases  where  Ph  is  very  small  limits  the  applications  of  these 
expressions.  However  by  using  conditional  probabilities  much  better  results  are  generated,  but 
the  correct  propagation  of  the  output  event  probabilities  is  a  major  factor  affecting  the  accu¬ 
racy  of  this  method.  To  derive  the  new  expressions,  the  expected  value  and  the  variance  are 
calculated  given  that  the  equivalent  edge  conducts.  Thus  the  following  expressions  are  pro¬ 
duced  from  (2.16)  and  (2. 17): 


Erfi  |  xW)  =  ,  Etx  I  x*0]E[y  I  y*0]  E[y  I  y*0]2V(x  I  x*0)  E[x  I  x*0]2  V(y  1  y*0)  2 

E[xlx*Oj  +  E[yly*0]  (E[x  I  x*0]  +  E[y  I  y*0])3  (E[xlx*0]  +  E[yly*0])3  ‘  * 


V(g  I  x*0,y*0)  =  [ - -WgQj - ]«  V(x  I  x*0)  +  [ - -[X-  — 1 - ]4 

v  1  r?r—  I .m  .  rr..  i  ..  iAi  J  i  r?r_.  i  ..  .  r-  r. .  t  J 


E[x  I  x*0]  +  E[y  I  y*0] 


E[x  1  x*0]  +  E[y  I  y*0] 


V(y)  (2.19) 


Where 


E[xlx*0]  =  -^L 

Phx 

(2.20) 

E[x2Ix*0]  = 

Phx 

(2.21) 

V(x  1  x*0)  =  E[x2 1  x*0]  -  E2[x  1  x^O] 

(2.22) 

This  method  produces  excellent  results,  but  unfortunately  it  is  dependent  on  the  accurate 
propagation  of  the  output  probabilities.  The  reason  that  the  second  method  works  better  than 
the  first  is  the  fact  that  the  approximation  for  the  first  method  is  hampered  by  the  presence  of 
a  large  probability  mass  at  zero.  Clearly  the  Gaussian  approximation  is  very  poor.  But  the 
conditional  distribution  will  follow  exactly  the  actual  distribution. 
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2.5  Calculation  of  Total  Current  and  Delay 

Having  devised  an  algorithm  for  graph  reduction  is  not  enough  for  calculating  the  current 
drawn  from  the  power  bus  by  the  circuit.  We  also  need  to  model  the  current  in  such  a  way 
that  the  information  that  we  get  from  the  graph  reduction  is  enough  to  fully  describe  the 
current  waveform  in  the  bus.  Following  the  implementation  in  [1]  the  current  in  the  bus  will 
be  described  by  two  triangular  pulses,  one  for  the  expected  value  and  one  for  the  variance. 
Therefore  a  peak  value  and  an  appropriate  time  span  must  be  calculated  in  both  these  cases. 
Also  the  output  event  must  have  some  delay  associated  with  it  for  proper  simulation  of  the 
circuit  In  this  chapter  the  method  used  in  II]  will  be  compared  with  a  new  method  that  has 
been  developed  in  this  thesis. 

The  basic  setup  that  is  going  to  be  used  is  shown  in  Fig.  2.6. 


CMOS  circuits  draw  current  from  the  bus  only  when  their  output  changes  state.  There¬ 
fore  the  calculations  for  the  current  must  be  done  when  an  input  is  switching.  The  method 
proposed  in  [1]  requires  the  calculation  of  the  probability  that  an  edge  is  conducting  given 
that  the  total  conductance  is  zero  for  every  edge  in  the  graph  that  is  being  reduced.  This 
requires  a  separate  graph  reduction  for  each  edge,  therefore  the  complexity  of  this  algorithm  is 
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0(n  ).  But  as  observed  in  [1]  this  is  not  very  accurate  because  we  cannot  be  certain  that 
imposing  the  condition  that  the  total  conductance  before  switching  occurs  is  zero  will  not 
affect  the  independence  of  the  inputs.  To  ensure  that  the  results  are  within  some  acceptable 
interval,  bounds  have  to  be  imposed  which  compromise  the  accuracy  of  the  method.  For 
completeness  we  present  the  expressions  for  the  current  statistics  as  they  are  derived  in  [1] 

E[i(t)]  -  E[ip(t)]  x  -  ~  +  E[i„(t)]  x  c  (2.23) 

CL  C 

E[i2(t)]  =  E[ip(t)]  x  (-r— 71—)2  +  E[in2(t)]  x  (■■  ■}  )2  (2.24) 

with  the  quantities  associated  with  E[ip(t)]  and  E[in(t)]  given  by  the  expressions 

E[ip]  =  Vdd  x  E[Gp(t+)  I  Gp(t~)  =  0]  x  P(Gp(D  =  0)  (2.25) 


E[ip2]  =  V2d  x  E[G2(t+)  I  Gp(r)  =  0]  x  P(Gp(D  =  0)  (2.26) 


These  are  the  results  for  the  p-block.  By  replacing  the  p  subscript  with  n  we  get  the 
results  for  the  n-part  as  well. 


The  new  algorithm  that  is  proposed  in  this  thesis  is  simpler  and  more  accurate  than  the 
one  just  described.  It  hias  been  assumed  that  the  inputs  are  independent  and  switch  instantane¬ 
ously.  Since  we  make  this  assumption,  then  the  probability  that  two  inputs  switch  at  the  same 
time  is  zero,  because  if  that  was  not  the  case  there  should  be  some  correlation  between  the 
signals.  Another  very  basic  observation  is  that  if  the  edge  that  switches  is  in  parallel  with  an 
edge  that  is  conducting,  then  switching  will  have  no  effect  whatsoever  to  the  output  as  illus¬ 
trated  in  Fig.  2.7. 
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This  edge 
is  ON 


This  edge 
switches 


Figure  2.7  Illustration  of  the  fundamentals  of  the  new  algorithm 

It  is  evident  that  the  edge  that  is  ON  effectively  shorts  the  edge  that  switches,  so  that  it 
has  no  effect  on  the  rest  of  the  circuit. 

Given  these  observations  we  can  state  that  for  the  output  to  switch,  one  of  the  edges 
must  switch  and  all  the  conducting  paths  that  are  created  must  go  through  the  edge  that  has 
switched.  Now  the  calculation  of  the  conductance  at  switching  becomes  very  simple.  Sup¬ 
pose  that  edge  k  is  switching.  Then  for  all  other  edges  i 

P(i  is  ON  I  k  is  switching)  =  P(i  is  ON)  =  Pw  (2.27) 

This  result  is  due  to  independence.  Since  at  switching  all  the  current  must  pass  through 
the  transistor  that  is  switching,  then  all  paths  that  bypass  that  transistor  should  not  conduct. 
To  ensure  this  during  elimination  we  use  the  following  algorithm.  In  the  description  of  the 
algorithm  we  split  the  edges  into  switching  and  non-switching.  A  switching  edge  is  defined 
as  the  edge  that  switches  and  every  equivalent  edge  that  contains  it.  It  is  evident  that  because 
of  independence  there  is  only  one  switching  edge  in  the  graph  that  is  being  reduced.  The 
only  edge  remaining  after  elimination  is  a  switching  edge.  Every  other  edge  is  a  non¬ 
switching  edge. 

1.  Assign  to  the  edge  that  switches  probability  1. 

2.  Proceed  to  the  elimination  of  the  graph  as  described  in  Section  2.3. 

2a.  Series  reduction  is  performed  as  usual. 
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2b.  In  the  parallel  reduction  we  distinguish  between  two  cases. 

-  If  two  non-switching  edges  are  in  parallel  then  elimination  proceeds  as 
usual. 


-  If  one  of  the  edges  is  switching  then  the  non-switching  edge  in  parallel 
with  it  is  assigned  Ph  =  0  and  the  probability  that  it  is  not  conducting  is 
stored.  Let  the  set  of  these  edges  be  denoted  by  a. 

3.  Calculate  E[I]  and  V(I)  from  the  following  expressions. 


and 


E[Ik]  =  Vdd  x  E[G(fromelimination)]  x  x  ]“|  Ph 

iea 


(2.28) 


E[Ik2I  =  V&  x  E[G2(ffomelimination)]  x  pihk  x  II  pii 

iea 


(2.29) 


The  philosophy  of  this  algorithm  is  that  calculating  the  equivalent  conductance  in  the 
way  that  was  just  described,  one  calculates  the  equivalent  conductance  given  that  the  circuit 
will  draw  some  current  from  the  bus.  Evidently  there  are  more  conducting  paths  than  the 
ones  calculated  by  this  method,  but  the  current  drawn  from  those  paths  at  the  time  of  switch¬ 
ing  is  zero  (internal  capacitances  are  considered  negligible).  So  this  method  calculates  the 
total  current  drawn  by  the  circuit.  The  expressions  just  mentioned  are  derived  very  simply 
through  the  assumption  of  independent  inputs,  because  the  result  of  the  elimination  is  the 
equivalent  conductance  given  that  an  edge  is  switching  and  the  edges  that  have  been  assigned 
probability  zero  are  not  conducting.  A  general  example  is  presented  in  the  following  Fig.  2.8. 
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Edges  that  are 
assigned 
probability  zero 
during  elimination 


Figure  2.8  Elimination  example 

The  algorithm  proposed  here  is  inherently  much  more  accurate  than  the  one  proposed  in 
[1]  and  is  also  0(n).  This  constitutes  a  great  improvement  over  the  previously  used  algorithm 
at  a  very  critical  step  towards  calculating  the  current 

Thus  far  the  peak  values  for  the  expected  value  and  variance  waveforms  have  been  cal¬ 
culated.  For  calculating  the  time  span  of  both  pulses  the  equations  derived  in  [1]  are 
sufficient  enough  for  the  proposed  approximation.  A  brief  description  follows. 

We  define: 

E[q]  =  +  V^CpPbk  (2.30) 

where  Pjho  and  Phlo  are  the  probabilities  that  the  output  will  go  from  "low"  to  "high"  and  vice 
versa.  Also 

E(k»  ■  E[y +  E[y  (2-3,) 

The  time  spans  for  the  expected  value  (x^)  and  the  variance  (xv)  are  given  by  the  follow¬ 
ing  simple  relations. 
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4  E[iq]  -  E[i]E[q] 

" "  3  (  V(i)  ’ 


(2.33) 


For  the  delay  of  the  output  event  a  weighted  average  of  the  delay  for  the  p-part  and  the 


delay  for  the  n-part  will  be  used.  Starting  from  the  deterministic  delay  model 


t,j  =  ln2  x 


Cp  +  Qn 


(2.34) 


and  using  simple  probability  theory  we  arrive  at  the  following  expressions  for  the  delay  of 
each  of  the  transitions. 


For  the  transition  from  low  to  high 


1  V(GP) 
*  (C’+C"> x  + 


(2.35) 


and  for  the  transition  from  high  to  low 


1  V(Gn) 

E(g  =  Mx(Vg,(— +  ~ 


(2.36) 


where  Gn  and  Gp  are  the  equivalent  conductances  of  the  n-part  and  the  p-part  when  they 
cause  switching.  For  the  total  delay  the  following  expression  can  be  used. 


E[td]  = 


E[t<in]  *  Eout,hl  E[tdp]  *  Eou^i}, 


Eout,hl  +  E0ut,lh 


(2.37) 


This  expression  combined  with  the  approach  that  was  used  in  its  derivation,  provides 
very  accurate  results,  which  are  also  easily  manipulated  (calculation  of  conditional  probabili¬ 
ties  and  other  related  quantities  is  greatly  simplified).  That  is  so  because  the  new  elimination 
process  breaks  down  the  graph  into  independent  edges  connected  in  series. 


In  [1]  some  different  expressions  had  been  proposed  but  the  straightforwardness  and  the 
simplicity  of  the  previous  expressions  in  conjunction  with  the  results,  introduce  this  method  as 
a  viable  candidate  for  the  calculation  of  the  delay.  The  above  expressions  are  simply  the 
deterministic  expressions  where  the  conductance  has  been  substituted  by  a  random  variable. 
According  to  [8]  this  is  the  best  second  order  approximation.  Further  tests  are  needed  to 
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determine  the  best  approach.  We  must  note  again  that  the  probabilities  used  in  the  previous 
expression  are  the  probabilities  given  that  the  output  is  switching  through  the  part  in  question. 

It  should  be  made  clear  that  this  analysis  is  for  one  channel  connected  partition  only  and 
not  for  the  entire  circuit.  Therefore  the  assumption  of  independence  is  valid  to  a  great  extent. 
As  was  demonstrated  in  [1]  independence  is  a  very  safe  approach.  In  the  case  of  correlated 
inputs  the  techniques  described  in  [1]  can  be  used,  but  the  stochastic  part  of  the  computation 
is  based  on  independence  and  this  is  what  has  been  presented  here.  As  for  the  assumption 
that  no  two  inputs  of  a  channel-connected  partition  can  switch  at  the  same  time,  it  is  valid  as 
long  as  the  inputs  are  independent.  But  this  can  easily  be  relaxed.  Let  us  assume  that  there 
are  several  signals  switching  at  the  same  time.  Because  of  independence  each  switching  input 
can  be  handled  separately  and  the  total  result  can  be  derived  by  superimposing  the  partial 
results.  It  is  very  interesting  to  point  out  that  this  has  already  been  done  in  [1]  with  excellent 
results.  In  any  case  this  is  not  an  essential  part  of  this  derivation  because  it  was  mainly  used 
to  avoid  the  cumbersome  notation  associated  with  two  edges  transitioning  at  exactly  the  same 
time,  which  was  of  no  practical  value  anyway. 

As  for  the  cases  where  two  current  pulses  are  overlapping  several  heuristics  can  be 
derived  in  order  to  approximate  the  real  current  waveforms. 

2.6  Results 

In  order  to  validate  the  methods  proposed  in  the  previous  chapters,  we  appied  the  new 
formulas  to  many  sample  circuits  and  compared  the  results  againsts  the  results  of  the  previous 
methods  and  the  results  of  a  SPICElike  simulation.  In  the  first  part  of  this  chapter  the  new 
reduction  scheme  is  compared  against  the  previous  one  and  in  the  latter  part  the  delay  calcula¬ 
tion  method  is  compared  with  the  exact  solution  given  by  the  SPICElike  simulation.  The 
gates  tested  include  both  simple  and  more  complex  structures. 


21 


First  examples  concerning  the  general  elimination  algorithm  will  be  displayed.  The 
results  will  be  compared  with  the  previously  existing  method  and  the  SPICElike  simulation. 
The  eliminated  graphs  range  from  three  transistor  to  multi-transistor  layouts.  Comparisons  for 
two  transistor  structures  (i.e.,  NAND  and  NOR  gate  configurations)  have  not  been  included 
because  of  their  simplicity. 

It  must  be  noted  that  "method  1"  is  the  method  in  which  no  signal  probabilities  are  used 
in  the  derivation  of  the  stochastic  properties  of  the  signals  (other  than  evaluating  the  output 
probability  and  the  characteristics  of  the  transistors  before  elimination).  The  comparison  is 
needed  to  calculate  the  effects  of  the  degradation  due  to  the  normal  approximation  without 
conditional  probabilities  as  explained  in  Section  2.3.  The  numbers  associated  with  the  arcs  on 
the  figures  correspond  to  the  probability  that  the  edge  conducts  and  its  maximum  conductance. 

"Method  2"  is  the  one  in  which  conditional  probabilities  are  taken  into  account  and,  as 
the  examples  will  demonstrate,  has  been  proven  the  most  accurate  among  the  three  approxi¬ 
mations. 

Example  1. 


0.4,40 


Figure  2.9  Example  1 
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E[g] 

V(g) 

CREST 

14.42 

39.70 

method  1 

14.97 

208.19 

method  2 

14.07 

297.23 

exact 

14.08 

297.29 

Table  2.1  Results  for  example  1 

This  simple  example  highlights  the  basic  notions  that  are  going  to  become  even  more 
apparent  in  later  examples.  First  of  all  method  2  is  much  more  accurate  than  CREST  in  com¬ 
puting  the  expected  value  and  the  variance  of  the  equivalent  edge.  This  is  something  con¬ 
sistent  throughout  the  comparison  of  the  two  techniques.  Secondly  method  1  is  quite  accurate 
but  not  as  accurate  as  the  other  two.  A  third  and  final  comment  is  that  CREST  underesti¬ 
mates  the  variance  by  a  great  amount  which  can  lead  in  turn  in  to  a  great  miscalculation  of 
the  current  characteristics.  As  to  the  particular  example,  method  2  has  shown  great 
effectiveness  in  calculating  this  particular  graph.  For  all  the  tested  values  for  the  probabilities 
as  well  as  the  conductances,  the  calculated  values  were  within  1  percent  of  the  values  pro¬ 
duced  by  the  SPICElike  "exact"  simulation. 
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Example  2. 


0.4 ,60  0.3 , 80 


0.7 ,40  0.5 , 100 

Figure  2.10  Example  2 


E[g] 

V(g) 

CREST 

21.71 

64.91 

method  1 

22.06 

254.28 

method  2 

20.26 

438.76 

exact 

20.41 

433.89 

Table  2.2  Example  2 

This  example  consists  of  four  transistors,  which  make  the  differences  more  apparent. 
The  error  for  method  2  is  much  better  than  that  of  CREST  both  for  the  expected  value  and 
the  variance.  Also  method  1  gives  reasonable  results  for  this  case. 
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Example  3. 


Figure  2.11  Example  3 


\ 

E[g) 

V(g) 

CREST 

4.54 

5.61 

method  1 

2.07 

199.70 

method  2 

4.15 

74.35 

exact 

4.17 

74.83 

Table  2.3  Example  3 

In  this  example  with  five  transistors  the  accuracy  difference  is  more  apparent.  Method  1 
fails  because  of  the  small  probabilities  involved  in  this  example  (i.e.,  the  0.1  and  0.2  conduct¬ 
ing  probabilities).  Apart  from  that  all  the  remarks  made  for  the  previous  examples  are  valid 

for  this  one  as  well. 

The  final  example  consists  of  thirteen  transistors  and  is  characteristic  of  the  behavior  of 


the  three  formulae. 
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Example  4. 


0.4,100 


E[g] 

V(g) 

CREST 

28.43 

38.50 

method  1 

29.46 

124.93 

method  2 

25.82 

332.74 

exact 

25.15 

349.50 

Table  2.4  Example  4 


The  final  conclusion  from  these  examples  is  that  method  2  is  consistently  better  than 
CREST  for  calculating  the  expected  value  and  the  variance.  The  added  precision  in  calculat¬ 
ing  the  expected  value  is  very  important  for  the  delay  and  current  estimation,  although  not  as 
spectacular  as  the  enhancement  in  the  calculation  of  the  variance,  which  is  also  very  impor¬ 
tant.  The  validity  of  method  1  has  been  demonstrated.  Although  it  is  not  optimal  it  does  not 
require  the  knowledge  of  the  edge  probabilities  once  the  stochastic  characteristics  of  each 
edge  have  been  established.  Finally  a  comprehensive  table  of  results  will  be  presented  to 
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better  judge  the  new  methods. 


#  of  edges 

m  2 

CREST 

m  1 

exact 

4 

5.29 

5.46 

3.55 

5.28 

5 

26.25 

30.76 

26.63 

26.66 

6 

2.41 

3.55 

diverges 

2.36 

8 

46.51 

49.70 

47.28 

46.62 

9 

54.98 

59.92 

54.59 

55.27 

10 

13.21 

14.82 

12.84 

13.27 

11 

10.97 

14.60 

7.19 

11.23 

12 

4.39 

4.88 

1.06 

4.40 

14 

55.24 

61.20 

55.17 

15 

47.00 

50.24 

56.90 

47.11 

17 

13.48 

15.13 

11.93 

11.53 

18 

93.90 

102.57 

89.89 

93.90 

19 

39.24 

45.83 

41.46 

39.38 

21 

55.68 

61.71 

48.25 

55.61 

23 

106.08 

115.80 

104.98 

109.10 

Table  2.5  Expected  value  examples 
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#  of  edges 

m2 

CREST 

m  1 

exact 

4 

273.0 

80.0 

666.2 

273.1 

5 

435.1 

254.7 

395.0 

393.3 

6 

49.8 

3.9 

diverges 

50.4 

8 

1282.8 

606.4 

1068.9 

1272.3 

9 

1329.5 

1232.6 

1325.6 

1335.9 

10 

243.1 

79.8 

268.4 

244.5 

11 

248.3 

118.3 

461.3 

255.5 

12 

78.6 

7.5 

150.3 

79.4 

14 

1522.5 

982.9 

1758.0 

1564.1 

15 

1308.3 

631.4 

890.9 

1297.0 

17 

249.4 

86.4 

300.5 

251.4 

18 

1914.6 

1401.9 

1838.4 

1994.3 

19 

353.5 

143.7 

192.6 

352.2 

21 

1543.7 

1009.1 

1859.1 

1585.3 

23 

2154.3 

1710.9 

2052.3 

2415.0 

Table  2.6  Variance  examples 

Finally,  we  give  a  simple  example  of  calculating  the  delay  and  the  equivalent  current. 
The  example  is  a  four  transistor  p-part  in  order  to  avoid  the  trivial  cases  of  two  or  three 
transistors  (for  which  the  aforementioned  quantities  can  be  calculated  exactly).  The  proba¬ 
bilistic  quantities  are  specified  so  that  only  the  p-part  will  conduct.  So  the  p-part  looks  like 
the  following  figure: 
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Figure  2.13  p-part  of  delay  and  current  example 
The  probabilistic  waveform  going  into  input  A  is  the  following: 


Figure  2.14  Input  waveform  for  example  in  Fig.  2.13 

From  the  analysis  the  following  results  have  been  produced  by  applying  (2.28)  to  (2.37): 

E[Ia]  =  1.557,  E[I i]  =  158.9,  E[delayA]  =  36.42  time  units 
The  actual  quantities  as  produced  by  exhaustive  SPICElike  analysis  are 

E[Ia]  =  1.567,  E[Ia]  =  160.1,  E [delay A]  =  35.96  time  units 
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3.  Layout  Extraction 

Layout  extraction  is  the  transformation  of  artwork  information  into  circuit  information. 
It  forms  a  very  important  link  between  the  physical  design  and  the  design  verification  process 
as  well  as  reliability  analysis.  In  particular,  computation  of  current  density  and  electromigra¬ 
tion  estimation  are  closely  linked  to  the  physical  process  and  the  layout  geometrical  informa¬ 
tion.  The  input  to  an  extractor  consists  of  the  mask  geometries  of  the  layout  together  with  a 
process  file.  The  artwork  is  typically  described  in  CIF  (Caltech  Intermedicate  Form),  an 
ASCII  format  where  each  layout  geometry  is  described  in  terms  of  its  shape  (usually  a  rectan¬ 
gle),  mask  layer,  and  coordinates.  The  output  is  typically  in  SPICE  format  which  can  be  read 
by  most  circuit  simulation  programs. 

In  the  previous  year,  we  developed  an  extractor  that  was  linked  to  the  OCT/VEM  design 
database  system  developed  at  the  University  of  California  at  Berkeley.  OCT  is  a  database 
system  capable  of  storing  the  design  for  an  entire  VLSI  chip,  including  the  schematic  artwork, 
and  the  extracted  netlist  of  a  design.  The  basic  unit  of  storage  is  a  cell,  which,  for  example 
may  contain  a  single  transistor,  a  logic  gate,  or  an  entire  functional  unit  (ALU,  CPU,  or 
RAM).  To  efficiently  represent  the  design  information,  cells  are  represented  hierarchically,  so 
that  a  cell  may  contain  instances  of  other  cells.  In  addition,  parameters  of  the  design  may  be 
annotated  to  the  structural  information,  such  as  the  dimensions  of  transistors,  or  parasitic  capa¬ 
citance  values  for  nets. 

VEM  is  the  graphical  editor  for  the  OCT  elements.  Circuit  schematics  of  artwork  may 
be  created  and  edited  with  VEM.  The  extractor  that  was  developed  during  the  first  year  of 
this  work  assumed  flattened  input  format  and  produced  flattened  SPICE  file.  However,  such  a 
flattened  representation  is  quite  inefficient  and  is  not  practical  for  very  large  designs  especially 
that  designs  are  specified  hierarchically.  During  this  year,  we  have  developed  a  new  extrac¬ 
tor,  based  on  the  previous  one.  The  new  extractor,  called  iCHARM,  preserves  and  exploits 
the  hierarchy  of  the  input  and  produces  a  hierarchical  circuit  file  as  well  as  the  power  bus 
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nets. 

* 

iCHARM  is  written  in  C,  runs  under  UNIX  and  has  the  following  features: 

CMOS  circuits: 

The  extractor  assumes  the  input  layout  is  in  CMOS  technology  with  any  number  of 
metal  interconnection  layers.  A  technology  file  is  read  to  determine  the  process  parame¬ 
ters  of  each  layer  and  the  connectivity  between  layers. 

Manhattan  geometries: 

iCHARM  supports  only  rectangles  as  mask  geometries;  therefore,  the  extractor  works 
only  for  manhattan-style  layouts;  i.e.,  layouts  where  the  edges  of  the  geometries  are 
parallel  to  the  x-  or  y-  axis. 

CIF  and  OCT  input  formats: 

The  input  layout  read  by  iCHARM  is  described  in  either  the  CIF  (Caltech  Intermediate 
Form)  format  or  in  the  OCT  database  format. 

Spice  output  format: 

The  extracted  circuit  description  is  output  in  the  SPICE  circuit  simulator  format.  If  the 
extraction  is  done  hierarchically,  each  cell  extracted  generates  a  .SUBCKT  card.  An 
option  to  output  the  circuit  in  the  OCT  format  also  is  envisioned. 

Scanline  algorithm: 

The  input  layout  is  represented  internally  in  iCHARM  in  terms  of  rectangles.  A  scanline 
algorithm  is  employed  to  carry  out  the  extraction  process.  This  is  described  in  detail  in 
this  thesis. 

The  program  may  be  run  in  several  modes  as  follows: 


UNIX  is  a  trademark  of  AT&T  Corp. 
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(1)  First,  the  user  has  the  option  of  running  the  extractor  in  flat  or  hierarchical  extraction 
mode.  If  there  is  a  high  degree  of  regularity  in  the  circuit  (if  the  circuit  has  just  a  few 
cells  but  perhaps  a  large  number  of  instances),  then  the  hierarchical  mode  should  run  fas¬ 
ter  and  use  less  run-time  memory  than  an  equivalent  run  in  the  flat  extraction  mode. 

(2)  The  user  also  has  a  choice  of  the  parasitics  of  the  interconnection  lines  that  may  be 
extracted.  The  default  mode  is  to  extract  each  net’s  substrate  capacitance  (the  capaci¬ 
tance  between  the  net  and  ground).  In  addition  to  the  capacitance,  the  user  may  extract 
the  resistance  through  each  interconnection  net. 

(3)  '  iCHARM  has  a  special  mode  to  extract  the  mask-level  rectangles  of  the  power  bus  nets 

(typically  Vdd  and  GND).  This  mode  does  additional  processing  of  the  power  bus  nets 
for  subsequent  analysis  by  other  programs  to  do  reliability  estimation. 

A  detailed  description  of  the  algorithms  used  in  iCHARM  together  with  the  implementation 
issues  are  included  in  a  M.S.  thesis  report  [4]. 

In  the  following,  we  describe  the  results  of  extracting  a  large  example,  the  2uchp  test 
chip,  specified  hierarchically  in  CIF. 

The  2uchp  design  is  a  10,000  transistor  CMOS  layout  that  was  designed  at  the  Rome  Air 
Development  Center  of  the  U.S.  Air  Force.  It  is  so  named  because  the  top-most  cell  is  named 
"2uchp."  The  design  was  originally  described  in  CIF,  but  it  was  converted  to  an  equivalent 
Oct  description  to  test  the  Oct  input  module  of  iCHARM. 

The  2uchp  chip  allowed  iCHARM  to  be  tested  with  realistic  layout  situations.  The  chip 
consists  of  74  different  cells.  The  design  hierarchy  of  the  chip  is  shown  in  Figure  3.1.  The 
major  cells  of  the  chip  are  a  ram  (rami),  rom  (rom2),  alu  ( alun2 ),  miscellaneous  logic 
(PRNG,  chp_qs,  prjeak),  and  the  padframe  (2ufrm).  Although  most  cells  also  contain  mask 
geometries,  in  Figure  3.1  only  the  subcells  are  shown  for  each  cell. 
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Figure  3.1  -  The 
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Since  the  majority  of  the  transistors  are  in  the  rami  cell,  it  is  the  most  difficult  cell  to 
extract  flat.  In  the  rami  cell  there  are  two  instances  of  the  memjblk  cell,  which  is  a  512-bit 
static  ram  array. 

The  2uchp  chip  contains  a  "hard"  extraction  case  where  cells  overlap:  the  2ufrm  cell  is  a 
padframe  cell  which  overlaps  the  other  major  cells  {rami,  rom2,  alunl,  PRNG,  chp_qs,  and 
pr_leak).  Another  major  cell,  the  alun2  cell,  has  overlapping  instances.  In  alun2,  the 
instance  of  alun2s  overlaps  the  other  instances  over  most  of  their  area.  Therefore  to  extract 
the  entire  chip  hierarchically,  it  is  important  to  avoid  flattening  the  areas  of  overlap. 

The  cells  of  the  chip  were  used  to  test  the  functionality  of  iCHARM.  All  of  the  cells 
were  put  through  the  flat  extraction,  parasitic  extraction,  power  bus  extraction,  and  hierarchi¬ 
cal  extraction  modes  of  the  extractor.  Of  most  interest  was  the  verification  of  the  hierarchical 
extraction  capability  and  the  performance  comparison  of  the  hierarchical  mode  versus  the  flat 
mode.  The  test  results  for  the  larger  cells  of  2uchp  are  shown  in  Table  3.1. 


Out  of  the  74  total  cells,  a  number  of  leaf  cells  are  simply  connector  or  feed-through 
cells  and  do  not  contain  any  transistors,  and  a  number  of  cells  are  simple  logic  gates  (nand, 
nor,  inv).  To  reduce  the  size  of  Table  3.1,  none  of  these  cells  are  reported.  The  second  and 
third  columns  of  Table  3.1  report  the  number  of  transistors  and  nets  for  each  cell.  Next,  the 
regularity  of  a  cell  is  defined  as 


regularity  = 


number  of  instances  of  the  cell 


number  of  unique  subcells 
For  instance,  a  cell  may  have  one  instance  each  of  five  subcells  and  four  instances  of  another 

subcell  which  results  in  a  regularity  factor  of  9/6  =  1.5.  The  number  of  overlap-rects  of  a  cell 

is  also  reported  in  Table  3.1. 


Cell  name 

Tran¬ 

sistors 

Table  3.1  -  iCHARM  Test  Results  with  2uchp  Cells 

Nets  Reg-  Overlap-  Flat  Hierarchical 

uiarity  rects  Time  Mem.  Time  Mem. 

col_sens 

16 

18 

50.0 

26 

0:00/0:00/0:03 

46.7 

0:00/0:01/0:04 

84.6 

ram_cell 

6 

8 

7.0 

0 

0:00/0:00/0:00 

9.5 

0:00/0:00/0:00 

12.8 

ram_cel4 

24 

19 

4.0 

0 

0:00/0:00/0:02 

29.9 

0:00/0:00/0:01 

29.1 

ram  nand 

8 

22 

20.5 

72 

0:00/0:00/0:01 

16.9 

0:00/0:00/0:01 

25.5 

ram  dec  na 

10 

7 

25.0 

47 

0:00/0:00/0:00 

15.1 

0:00/0:00/0:01 

24.5 

ram_rowa 

28 

39 

2.0 

57 

0:00/0:00/0:02 

46.4 

0:00/0:01/0:03 

69.6 

mem_blk 

3072 

1122 

58.7 

1 

0:00/0:07/4:47 

3022.8 

0:00/0:00/0:18 

945.7 

ram_rowd 

448 

145 

18.5 

18 

0:01/0:02/0:39 

502.0 

0:01/0:01/0:11 

393.3 

alun2s 

0 

4 

3.0 

301 

0:00/0:00/0:00 

1.8 

0:00/0:00/0:00 

2.6 

tri  8 

14 

16 

40.0 

1 

0:00/0:00/0:01 

28.3 

0:00/0:00/0:02 

57.0 

rom_dec  el 

4 

8 

9.0 

10 

0:00/0:00/0:00 

4.4 

0:00/0:00/0:00 

6.0 

inv  add 

6 

4 

11.5 

16 

0:00/0:00/0:00 

7.9 
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14.1 
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24 

26 

70.0 

51 
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27.1 
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39.6 
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2 

7 

2.0 

12 
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2.4 
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3.6 
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305 

80.0 

370 
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2 

6 
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6 
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3.0 
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55 

17.7 

49 
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106.3 
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134.9 
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23 

5 
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4 
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96.4 
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149.7 
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34 

10 

215.0 

1 

0:01/0:01/0:04 

124.7 

0:01/0:03/0:09 

179.4 
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16 

5 
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1 
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0:00/0:01/0:04 
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8 
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5 
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8.9 

0:00/0:00/0:00 

15.1 
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9.0 

1379 
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2429.7 

0:04/0:09/3:32 

6769.8 
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18 

27 

1.3 

58 

0:00/0:00/0:02 

33.0 

0:00/0:00/0:02 

54.3 
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50 

53 

1.5 

69 
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108.1 
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29 

?? 

O.M. 
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25.9 
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?? 
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Input  format:  Oct  format 
Hardware:  Vax  3500 

Time  (user)  reported  in  units  of  minutes:seconds,  see  text  for  format 
O.M.  =  Out  of  Memory  (the  process  aborted) 


Parasitics  extracted:  Capacitance  only 
Memory  reported  in  units  of  kbytes 
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The  performance  results  for  flat  and  hierarchical  extractions  are  reported  in  the  form: 
input/pre-process/final.  All  reported  times  are  the  user  times  taken  by  the  process  after  the 
given  step  has  completed.  The  input  time  refers  to  the  time  taken  to  read  the  layout  in. 
Some  preprocessing  is  done  before  the  actual  extraction;  the  time  after  the  preprocessing  step 
is  reported.  In  flat  extraction,  the  preprocessing  step  is  the  flattening  of  the  layout,  and  in 
hierarchical  extraction,  the  overlap  analysis  is  the  preprocessing  step.  The  final  time  is  the 
total  time  taken  by  the  process.  In  other  words,  all  times  are  reported  relative  to  the  start  time 
of  zero,  therefore,  "0:04/0:12/2:13"  means  the  process  took  4  seconds  of  input  time,  then  8 
seconds  for  preprocessing,  followed  by  2:01  minutes  to  do  the  extraction. 

Since  the  memory  allocation  in  iCHARM  was  handled  manually,  it  was  possible  to  moni¬ 
tor  the  amount  of  memory  that  was  allocated  for  each  run  of  the  program.  The  maximum 
amount  of  memory  that  was  allocated  during  the  process  is  reported  in  Table  3.1. 

In  looking  over  the  test  results  in  Table  3.1,  some  interpretations  and  conclusions  can  be 
formed  from  the  data. 

First  of  all,  the  overhead  in  the  preprocessing  step  to  handle  overlap  is  minimal.  In  most 
cases  the  time  to  process  the  overlap  is  less  than  the  time  to  flatten  the  entire  cell.  A  good 
example  of  this  is  the  2ufrm  cell,  where  it  took  12  seconds  to  flatten  but  only  9  seconds  to 
process  the  overlap.  Most  cells  in  Table  3.1  take  almost  negligible  time  to  process  the  over¬ 
lap:  from  "zero"  (less  than  1  second)  to  2  seconds. 

To  compare  the  times  for  hierarchical  extraction  versus  flat  extraction,  it  makes  sense  to 
discuss  only  the  cells  that  have  runtimes  greater  than  10  seconds.  Measurement  errors  for  the 
runtimes  less  than  10  seconds  are  a  large  percentage  of  the  total  runtime,  so  the  shorter  results 
cannot  be  trusted  for  comparisons.  Of  the  results  in  Table  3.1,  there  are  18  cells  where  the 
total  flat  extraction  time  took  longer  than  10  seconds.  These  cells  can  be  classified  into  three 


groups. 


The  best  results  show  the  advantages  of  hierarchical  analysis:  the  hierarchical  extraction 
takes  less  time  and  less  memory.  This  occurs  for  the  largest  number  of  cells,  1 1  cells  in  all: 
[memjblk,  ram_rowd ,  col_blk_r ,  coljblkj,  alur3,  alur2 ,  deccol,  PRNGb,  PRNGa,  rami , 
2uchp }.  The  extraction  has  a  few  overlap-rects,  which  means  there  is  little  or  no  overlap  for 
the  cell.  In  fact,  the  memjblk  cell  has  virtually  no  overlap;  therefore,  the  speed-up  is  almost 
16-fold,  and  one-third  of  the  memory  is  used.  For  the  rami  and  2uchp  cells,  hierarchical 
extraction  produces  an  extracted  result  where  the  flat  analysis  failed  to  produce  any  result. 

When  overlap  is  present,  the  results  are  less  impressive;  five  cells  still  have  faster 
hierarchical  extraction  times  but  more  memory  is  actually  used:  [rom_ar2,  data,  ram_cold, 
rom2,  PRNG).  All  of  the  cells  end  up  having  significant  lists  of  overlap-rects.  More  memory 
must  be  used  to  calculate  and  create  the  overlap-rects.  The  extraction  process  itself  is  slowed 
since  during  extraction  each  rect  for  the  cell  must  be  checked  to  see  if  it  touches  an  overlap- 
rect.  Nevertheless,  the  savings  due  to  hierarchy  are  great  enough  so  that  the  runtime  is  still 
faster. 

Finally  there  are  two  cells  for  which  hierarchical  analysis  is  worse  in  terms  of  both  speed 
and  memory  usage:  [2ufrm,  alun2 } .  As  mentioned  before  the  padframe  is  the  "hard”  test-case 
for  hierarchical  extractors,  and  the  alun2  cell  contains  a  significant  amount  of  overlap.  How¬ 
ever,  it  does  appear  that  for  the  18  test-cases,  most  of  the  cells  benefit  from  being  extracted 
hierarchically  with  the  methods  employed  by  iCHARM. 
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4.  Estimation  of  Current  Density  in  the  Bus 

4.1  Introduction 

One  of  the  primary  aims  of  this  work  is  the  estimation  of  the  current  density  in  any  sec¬ 
tion  of  the  bus  (power  or  ground)  for  electromigration  estimation  and  for  re-sizing  the  bus  for 
electiomigration  prevention,  or  at  least  increasing  the  MTF  due  to  electromigration. 

During  the  first  year  of  this  work,  we  developed  a  computer  program  for  extracting  the 
RC  model  of  the  bus  and  for  computing  the  average  current  densities  in  the  bus.  The  bus 
extractor  is  based  on  decomposing  the  bus  geometry  into  geometrical  primitives,  such  as 
straight,  L,  and  T  primitives.  Each  primitive  is  pre-characterized  using  finite-element- 
methods.  During  extraction  of  the  bus  the  RC  model  of  each  extracted  primitive  is  used.  A 
linear  solver  is  then  used  to  solve  for  the  average  current  density  in  any  section  of  the  RC  bus 
nodal,  given  the  average  current  density  at  each  contact  node  as  obtained  from  the  probabilis¬ 
tic  simulator. 

During  the  second  phase  of  this  work,  we  improved  the  analysis  of  the  bus  in  three 
directions.  First,  we  completed  the  modeling  of  the  contacts.  Strips  of  metal  containing  con¬ 
tacts  were  analyzed  using  the  finite  element  method  for  various  combinations  of  contact  sizes 
and  metal  widths,  and  the  result  store  in  the  set  of  primitives. 

Second,  we  have  shown  that  MTF  due  to  electromigration  depends  not  only  on  the  aver¬ 
age  current  density,  but  also  on  the  variance  of  the  current  density  waveform  [2].  For  exam¬ 
ple,  two  density  waveforms  having  the  same  average,  but  different  variance  may  lead  to 
different  MTF  estimation.  This  fact  has  been  observed  experimentally  and  reported  in  litera¬ 
ture. 

The  third  issue  considered  is  the  computation  of  the  variance  current  density  waveform 
in  the  bus,  in  addition  to  the  average  current  density  waveform.  Recall  that  the  variance 
current  waveforms  drawn  by  the  logic  gates  at  the  bus  contacts  are  computed  by  the  proba- 
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bilistic  simulation,  in  addition  to  the  average  current  waveform.  Earlier,  we  have  only  con¬ 
centrated  on  the  computation  of  the  average  current  density  [5]  and  the  calculation  of  the  vari¬ 
ance  was  not  done.  In  this  work  we  have  developed  algorithms  and  software  for  efficient 
computation  of  both  the  average  and  the  variance  current  densities  in  the  bus. 

The  problem  of  finding  average  and  variance  of  branch  currents  in  a  general  RC  network 
was  first  published  by  Najm  et.  al.  in  [2].  The  method  involves  finding  the  infinite  impulse 
response  functions  for  each  of  the  branch  currents  of  the  network  due  to  unit  current  impulses 
applied  at  contact  points,  one  at  a  time.  To  calculate  the  variance  of  branch  currents,  one 
needs  to  carry  out  a  convolution  operation  involving  infinite  terms.  Such  an  analysis  is 
extremely  complex  and  perhaps  not  necessary  when  the  capacitive  effects  of  buses  are  almost 
negligible  [9].  Even  for  simple  resistive  networks,  the  above  method  involves  inverting  the 
node  admittance  matrix,  which  is  very  expensive  both  from  computation  and  storage  points  of 
view.  In  this  paper,  we  have  simplified  the  above  analysis  for  resistive  networks  and  have 
made  it  very  efficient.  The  sparsity  of  the  node  admittance  matrix  is  exploited  to  reduce  CPU 
time  and  memory  space  requirements. 

We  assume  that  there  are  n  nodes  and  b  branches  in  the  resistive  network  that  represents 
the  bus.  Let  lb  denote  the  vector  of  branch  currents.  Is  the  vector  of  currents  entering  the 
nodes  (contact  point  currents)  and  Vn  the  vector  of  node  voltages.  We  denote  j*  element  of  a 
vector  V  by  Vj  and  the  ij*  element  of  a  matrix  A  by  Ay.  The  i*  column  vector  of  a  matrix 
A  will  be  denoted  by  [A]j.  The  transpose  of  a  matrix  A  will  be  denoted  by  AT.  We  assume 
that  the  contact  point  current  waveforms  (both  average  E[Is]  and  variance  Var[Is])  have 
already  been  sampled  in  time  and  there  are  m  time  points  for  which  the  network  needs  to  be 
analyzed.  Furthen  ore,  as  justified  in  [1],  we  assume  that  various  contact  point  currents  are 
mutually  independent  at  a  given  time  point. 
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4.2  Networks  containing  no  loops  (TREE): 

For  a  network  with  no  loops,  current  in  a  branch  (say  k*  one)  can  be  expressed  as  a 
linear  combination  of  the  leaf  (contact)  currents: 

n 

Ibk  =  £  hjk  ISj  (4.1) 

j=l 

where  h^  e  0,1.  h^  =  1  if  node  j  is  a  descendent  of  branch  k  and  =  0  otherwise.  By  taking 
expectation,  we  get  the  average  value  of  the  branch  current: 

EM  =  X  hjk  E[Isj].  (4.2) 

J=i 

Since  various  contact  currents  are  assumed  to  be  mutually  independent,  we  can  write: 

Var[IbjJ  =  £  hjk  Var[ISj]  =  £  Var[ISj],  (4.3) 

j=i  j=i 

Therefore,  by  analyzing  the  tree  is  in  a  bottom  up  fashion  (from  contact  points  up  to  the  pad), 
we  obtain  the  average  and  the  variance  of  currents  in  various  branches.  On  a  serial  computer, 
this  algorithm  would  take  O(mn)  time  in  the  worst  case. 


4.3  Networks  containing  loops  (GENERAL): 

For  a  general  resistive  network,  to  calculate  branch  currents,  we  need  to  calculate  node 
voltages  first.  The  following  equation  holds  between  Vn  and  Is,  where  Yn  is  the  node  admit¬ 
tance  matrix  of  the  network: 

Yn  Vn  =  Is.  (4.4) 

By  taking  expectation: 

Yn  E[Vn]  =  E[Is].  (4.5) 

From  the  LU  factors  of  the  Yn  matrix,  this  equation  can  be  solved  to  obtain  the  average 
values  of  node  voltages.  If  we  assume  that  i*  resistor  (conductance  =  y4)  is  connected 
between  nodes  1  and  k  in  the  network,  then 
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E[Ibj]  =  (E[VnJ  -  E[Vnk])  x  yr  (4.6) 

Thus,  we  can  calculate  the  average  values  of  currents  in  all  the  branches.  The  calculation  of 
variance  is  more  complex.  We  define  the  following  covariance  matrices: 

KIb  4  E[(Ib  -  E[Ib])(Ib  -  E[Ib])T  (4.7a) 

KIs  4  E[(Is  -  E[Is])(Is  -  E[Is])t]  (4.7b) 

KVn  4  E[(Vn  -  E[Vn])(Vn  -  E[Vn])T]  (4.7c) 


The  diagonal  entries  of  KIb  are  the  required  variance  entries  of  interest.  As  before,  if  the  i* 
resistor  (conductance  =  y,)  is  connected  between  nodes  1  and  k,  then  variance  of  current  in 
the  resistor  is: 

Klb^  =  (KVnu  +  KVn,*  -  2KVnlk)  x  yi2.  (4.8) 

Thus,  to  calculate  the  diagonal  entries  of  KIb,  only  (n+b)  elements  of  KVn  matrix  are  needed. 
Elements  needed  are  the  n  diagonal  elements  and  b  off  diagonal  elements  (KVn*)  correspond¬ 
ing  to  each  resistor  between  nodes  1  and  k.  (Note  KVnlk  =  0  if  1  or  k  is  the  reference  node.) 

To  find  the  required  KVn  entries,  we  proceed  as  follows.  Eq.  (4.4)  can  be  written  as 
Vn  =  Rn  Is.  Rn  is  the  resistance  matrix,  which  is  inverse  of  Yn.  By  substituting  this  in  Eq. 
(4.7c),  we  get  the  following  relationship  between  KVn  and  KIs: 

KVn  =  Rn  KIs  Rn.  (4.9) 

Note  that  Rn  is  a  symmetric  matrix.  KIs  is  a  diagonal  matrix  because  of  the  assumption  that 
various  contact  currents  are  mutually  independent.  All  of  the  diagonal  entries  of  KIs  are  posi¬ 
tive,  i*  entry  being  Var(lsi)  (=  els2  say).  Thus,  KIs  is  a  positive  definite  matrix  and  there¬ 
fore  it  has  a  positive  definite  diagonal  square  root  (call  it  als): 

KIs  =  als  als  =  (als)2-  (4.10) 

The  i*  diagonal  entry  of  als  is  alSj.  This  equation  can  be  used  to  express  KVn  matrix  as  fol¬ 
lows  (since  Rn  is  symmetric): 
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K Vn  =  (Rn  aIs)(Rn  oIs)T-  (4. 1 1 ) 

Let  aVn  =  Rn  als.  Note  that  aVn  is  not  symmetric.  Equation  (4.11)  becomes: 

KVn  =  (cVn)(oVn)T-  (4.12) 

If  we  know  aVn,  we  can  calculate  scriptl.k*  entry  of  KVn  by  the  equation: 

n 

KVnik  =  Z  oVnij  x  oVnkj-  (4.13) 

j=i 

Since  Yn  is  the  inverse  matrix  of  Rn,  we  can  rewrite  oVn  =  Rn  als  as: 

Yn  aVn  =  als.  (4.14) 

For  each  of  the  n  corresponding  columns  of  aVn  and  als  matrices,  we  get  n  matrix  equations 
of  the  form:  (i  =  1,  2,  —  n) 

>  r  ■s 

aVnH  0 

aVn2i  : 

Yn  ;  =  alSj  or  Yn[aVn]j  =  [als^.  (4.15) 

aVn2i  : 

0 

One  can  solve  these  n  equations  and  obtain  the  required  KVn  entries  from  Eq.  (4.13).  How¬ 
ever,  the  node  admittance  matrix  of  a  resistive  network  is  usually  very  sparse  and  therefore, 
even  for  large  values  of  n,  it  can  be  stored  very  compactly  on  a  computer.  On  the  other  hand, 
aVn  matrix  is  usually  a  full  matrix.  Therefore,  for  large  values  of  n,  storage  would  be  a 
problem.  However,  from  Eq.  (4.13),  we  notice  that  to  calculate  KVnlk,  at  any  time,  we  do  not 
need  to  multiply  two  entries  of  aVn  matrix  from  two  different  columns.  Therefore,  we  do  not 
need  to  compute  and  store  the  entire  aVn  matrix  at  any  time.  As  we  calculate  columns  of 
oVn,  we  perform  the  multiplications  needed  in  Eq.  (4.13)  for  all  the  required  KVn^  entries 
and  then  discard  the  column  vector.  If  we  assume  that  the  solution  phase  (backward  and  for¬ 
ward  substitutions)  requires  0(n“)  multiplications  then  to  calculate  average  and  variance 
current  waveforms,  we  need  0(n1+a  +  mn2  +  mn(n+b))  multiplications.  Computational 
speedup  can  be  obtained  by  using  parallel  processing  since  the  above  algorithm  is  quite  suit- 
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able  for  implementation  on  a  medium  to  coarse  grain  parallel  machine. 

4.3  Experimental  Results 

In  this  section,  we  demonstrate  the  effectiveness  of  the  various  algorithms  described  ni 
the  previous  section  on  several  examples.  All  of  the  experiments  reported  in  this  section  were 
carried  out  on  a  VAX  3500  computer.  Since  the  main  emphasis  of  this  part  of  the  work  has 
been  to  come  up  with  efficient  algorithms  for  the  calculation  of  average  and  variance  of 
branch  currents,  we  have  randomly  generated  several  resistive  networks  and  have  assigned 
random  values  for  average  and  variance  of  contact  currents.  In  Table  I,  we  report  results  of 
the  experiments  for  TREE  networks.  The  pre-processing  phase  consists  of  generation  of  the 
network  and  setting  up  of  various  data  structures.  Table  II  shows  similar  results  for  GEN¬ 
ERAL  networks.  For  the  LU  decomposition  and  the  solution  phase,  we  have  used  the  sparse 
linear  equation  solver  package  [10]. 

The  details  of  the  derivation  of  the  above  solution  technique,  with  examples,  are  given  in 


Ref.  [6]. 
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V.  Summary 

In  this  report,  we  have  summarized  our  accomplishments  over  the  past  year  on  the  prob¬ 
lem  of  computing  current  density  statistics  in  CMOS  circuits  for  reliability  estimation.  The 
emphasis  has  been  on  the  use  of  the  average  value  and  variance  of  the  current  density  to  esti¬ 
mate  MTF  due  to  electromigration  in  the  bus.  The  work  involves  three  subtasks:  (1)  extrac¬ 
tion  from  layout;  (2)  probabilistic  simulation  of  transistor  circuit  to  compute  the  average  value 
and  variance  of  current  drawn  by  the  circuit  at  contact  points;  and  (3)  the  computation  of  the 
current  density  in  the  bus  and  estimation  of  the  MTF.  An  overview  of  the  system  that  we 
developed  is  given  in  [7]. 

In  the  work  on  extraction,  we  have  developed  a  hierarchical  extractor  based  on  our  previ¬ 
ous  extractor.  The  new  extractor  preserves  the  hierarchy  inherent  in  the  design,  thus  allowing 
hierarchical  simulation  and  saving  storage.  In  probabilistic  simulation,  we  have  extended  and 
enhanced  the  accuracy  of  our  initial  derivation  and  started  the  development  of  a  hierarchical 
probabilistic  simulator  that  exploits  the  hierarchy  in  the  design  and  the  data  structure.  For  the 
computation  of  the  current  density  in  the  bus,  we  have  developed  circuit  reduction  and  sparse 
matrix  techniques  to  compute  the  average  value  as  well  as  the  variance  current  density.  The 
RC  circuit  model  of  the  bus,  although  linear,  is  very  large  that  could  contain  more  than  one 
hundred  thousand  nodes.  Special  numerical  techniques  are  then  necessary  to  handle  the  com¬ 
plexity  of  the  computation. 

So  far,  we  have  concentrated  our  efforts  on  estimating  the  current  density  for  estimating 
MTF  due  to  electromigration.  However,  we  believe  that  the  probabilistic  approach  is  applica¬ 
ble  to  other  reliability  studies  as  well.  We  have  thus  initiated  two  projects  in  that  direction. 
One  is  the  estimation  of  hot-carrier  effects  in  transistors.  Hot-carrier  degradation  is  also  a 
long-term  reliability  issue  that  is  cumulative  in  nature,  depending  on  how  long  a  transistor 
operates  in  saturation  under  all  possible  primary  input  conditions,  and  thus  depends  on  its 
location  in  the  circuit.  Exhaustive  simulation  would  be  exorbitantly  expensive  to  estimate 


45 


hot-carrier  effects.  We  are  currently  developing  probabilistic  techniques  to  estimate  these 
effects  with  reasonable  computational  effort  We  are  also  developing  techniques  to  estimate 
voltage  drop  in  the  busses. 

Voltage  drop  is  also  becoming  of  concern  with  the  decrease  in  feature  size.  It  will 
become  more  critical  when  the  supply  voltage  is  reduced  and  when  BICMOS  circuits  are  used 
in  the  design  since  BICMOS  circuits  draw  a  larger  current  during  switching  that  purely 
CMOS  circuits.  In  contrast  to  electromigration  and  hot-carrier  effects,  voltage  drop  is  a 
worst-case  phenomenon,  and  our  aim  is  to  develop  techniques  to  sort  out  these  worst-case 
situations.  Our  ultimate  aim  is  to  optimize  a  given  design,  or  give  design  guidelines,  to  keep 
electromigration  and  voltage  drop  within  acceptable  bounds  everywhere  in  the  bus,  as  well  as 
keep  degradation  with  transistors  caused  by  hot-carrier  effects  within  operating  limits,  for  all 
possible  inputs. 
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